***Table 1***
clear
use "StateData.dta"
order sip_lead7 sip_lead6 sip_lead5 sip_lead4 sip_lead3 sip_lead2 sip_lead1 sip sip_lag* newcases_lag* newdeaths_lag* mobility_lag* unemployment_lag*
reghdfe newcases sip newcases_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using Table1, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
foreach dv in newdeaths mobility unemployment {
reghdfe `dv' sip `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using Table1, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
}


***Figure 1***
clear
use "StateData.dta"
egen maxsip = max(sip), by(state)
g never = maxsip < 1
g sipday = day if sip == 1
egen maxsipday = max(sipday), by(state)
g stayon = maxsipday == 22066
g onoff = never + stayon == 0
drop max*
*recode day as Days after Feb 24
replace day = day - 21969
foreach i in AL AR AZ CA CO CT DC DE FL GA HI IA ID IL IN KS KY LA MA MD ME MI MN ///
MO MS MT NC ND NE NH NJ NM NV NY OH OK OR PA RI SC SD TN TX UT VA VT WA WI WV WY {
preserve
keep if state == "`i'"
lpoly mobility day, bw(7) nograph gen(k_mobility) at(day)
sum never
if r(mean) == 1 {
	sort day
	graph twoway (line k_mobility day) (scatter mobility day)
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush Mobility
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule -.75 .25 .25 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore
	}
sum stayon
if r(mean) == 1 {	
	sum day if sip == 1
	local start = r(min)
	sort day
	graph twoway (line k_mobility day) (scatter mobility day) , xline(`start')
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush Mobility
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule -.75 .25 .25 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore	
	}
sum onoff
if r(mean) == 1 {
	sum day if sip == 1
	local start = r(min)
	local end = r(max)
	sort day
	graph twoway (line k_mobility day) (scatter mobility day) , xline(`start') xline(`end')
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush Mobility
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule -.75 .25 .25 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
	gr_edit .plotregion1._xylines[2].style.editstyle linestyle(color(black) pattern(dot)) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore
	}	
}
graph combine "AL.gph" "AR.gph" "AZ.gph" "CA.gph" "CO.gph" "CT.gph" "DC.gph" ///
"DE.gph" "FL.gph" "GA.gph" "HI.gph" "IA.gph" "ID.gph" "IL.gph" "IN.gph" "KS.gph" ///
"KY.gph" "LA.gph" "MA.gph" "MD.gph" "ME.gph" "MI.gph" "MN.gph" "MO.gph" ///
"MS.gph" "MT.gph" "NC.gph" "ND.gph" "NE.gph" "NH.gph" "NJ.gph" "NM.gph" ///
"NV.gph" "NY.gph" "OH.gph" "OK.gph" "OR.gph" "PA.gph" "RI.gph" "SC.gph" ///
"SD.gph" "TN.gph" "TX.gph" "UT.gph" "VA.gph" "VT.gph" "WA.gph" "WI.gph" ///
"WV.gph" "WY.gph", col(5)
gr_edit .style.editstyle declared_ysize(7) editcopy
gr_edit .style.editstyle margin(zero) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
graph export "Figure1.pdf", replace as(pdf)


***Figure 2***
*Newcases
clear
use "StateData.dta"
reghdfe newcases sip sip_lag* newcases_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on New Cases per Million Residents
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 14 2 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph save "LaggedEffects_NewCases.gph", replace
*Deaths
clear
use "StateData.dta"
reghdfe newdeaths sip sip_lag* newdeaths_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Deaths per Million Residents
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 14 2 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph save "LaggedEffects_Deaths.gph", replace
*Mobility
clear
use "StateData.dta"
reghdfe mobility sip sip_lag* mobility_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Mobility
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 14 2 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph save "LaggedEffects_Mobility.gph", replace
*Unemployment
clear
use "StateData.dta"
reghdfe unemployment sip sip_lag* unemployment_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Insured Unemployment Rate
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 14 2 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph save "LaggedEffects_Unemployment.gph", replace
*Combine into one figure
graph combine "LaggedEffects_NewCases.gph" "LaggedEffects_Deaths.gph" "LaggedEffects_Mobility.gph" "LaggedEffects_Unemployment.gph", col(2)
gr_edit .plotregion1.graph1.yaxis1.title.text = {}
gr_edit .plotregion1.graph1.yaxis1.title.text.Arrpush Effect of SIP
gr_edit .plotregion1.graph3.yaxis1.title.text = {}
gr_edit .plotregion1.graph3.yaxis1.title.text.Arrpush Effect of SIP
gr_edit .plotregion1.graph1.xaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph2.xaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph2.yaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph4.yaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph1.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph1.title.text = {}
gr_edit .plotregion1.graph1.title.text.Arrpush Cases
gr_edit .plotregion1.graph2.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph2.title.text = {}
gr_edit .plotregion1.graph2.title.text.Arrpush Deaths
gr_edit .plotregion1.graph3.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph3.title.text = {}
gr_edit .plotregion1.graph3.title.text.Arrpush Mobility
gr_edit .plotregion1.graph4.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph4.title.text = {}
gr_edit .plotregion1.graph4.title.text.Arrpush Unemployment
gr_edit .style.editstyle margin(zero) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
graph export "Figure2.pdf", replace as(pdf)


***Tables A1-A4***
clear
use "StateData.dta"
order sip_lead7 sip_lead6 sip_lead5 sip_lead4 sip_lead3 sip_lead2 sip_lead1 sip sip_lag* newcases_lag* newdeaths_lag* mobility_lag* unemployment_lag*
foreach dv in newcases newdeaths mobility unemployment {
reghdfe `dv' sip, a(stateid day) cluster(stateid)
qui:outreg2 using table_`dv', replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip_lead* sip , a(stateid day) cluster(stateid)
qui:outreg2 using table_`dv', append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip sip_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using table_`dv', append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using table_`dv', append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip_lead* sip sip_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using table_`dv', append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip_lead* sip `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using table_`dv', append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip sip_lag* `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using table_`dv', append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip_lead* sip sip_lag* `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using table_`dv', append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
}


***Table A5***
clear
use "StateData.dta"
sort state day
g change_sip = sip - sip[_n-1] if state == state[_n-1] & day == day[_n-1] + 1
table day, content(mean change_sip min change_sip max change_sip)
*first pullback of sip is on April 27, only analyze days before that
drop if day >= 22032
reghdfe newcases sip newcases_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA5, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe newdeaths sip newdeaths_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA5, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe mobility sip mobility_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA5, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe unemployment sip unemployment_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA5, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)


***Table A6***
clear
use "StateData.dta"
sort state day
g change_sip = sip - sip[_n-1] if state == state[_n-1] & day == day[_n-1] + 1
table day, content(mean change_sip min change_sip max change_sip)
*last expansion of sip is on April 7, only analyze after that
drop if day <= 22012
reghdfe newcases sip newcases_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA6, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe newdeaths sip newdeaths_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA6, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe mobility sip mobility_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA6, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe unemployment sip unemployment_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA6, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)


***Table A7***
clear
use "StateData.dta"
egen maxpop = max(pop), by(state)
sum maxpop, d
g smallstate = maxpop < r(p50)
tab state smallstate
keep if smallstate == 1
reghdfe newcases sip newcases_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA7, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe newdeaths sip newdeaths_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA7, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe mobility sip mobility_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA7, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe unemployment sip unemployment_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA7, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)


***Table A8***
clear
use "StateData.dta"
g sipday = day if sip == 1
egen minsipday = min(sipday), by(state)
replace minsipday = 22012 if minsipday == .
sum minsipday, d
g laterstate = minsipday > r(p50)
tab state laterstate
keep if laterstate == 1
reghdfe newcases sip newcases_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA8, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag14)
reghdfe newdeaths sip newdeaths_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA8, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag14)
reghdfe mobility sip mobility_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA8, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag14)
reghdfe unemployment sip unemployment_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA8, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag14)


***Table A9***
clear
use "StateData.dta"
reghdfe newcases sip neighbor newcases_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA9, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe newdeaths sip neighbor newdeaths_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA9, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe mobility sip neighbor mobility_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA9, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe unemployment sip neighbor unemployment_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA9, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)


***Table A10***
clear
use "CountyData.dta"
reghdfe newcases sip newcases_lag*, a(fips day) cluster(stateid)
qui:outreg2 using TableA10, replace dec(3)  alpha(0.01, 0.05) sortvar(sip)
foreach dv in newdeaths mobility {
reghdfe `dv' sip `dv'_lag*, a(fips day) cluster(stateid)
qui:outreg2 using TableA10, append dec(3)  alpha(0.01, 0.05) sortvar(sip)
}


***Table A11***
clear
use "StateData.dta"
reghdfe newcases sip newcases_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA11, replace dec(3)  alpha(0.01, 0.05) sortvar(sip)
reghdfe newdeaths sip newdeaths_lag*, a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA11, append dec(3)  alpha(0.01, 0.05) sortvar(sip)


***Table A12***
clear
use "StateData.dta"
order sip_lessthan2weeks sip_2weeks
reghdfe newcases sip_lessthan2weeks sip_2weeks newcases_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA12, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lessthan2weeks sip_2weeks)
foreach dv in newdeaths mobility unemployment {
reghdfe `dv' sip_lessthan2weeks  sip_2weeks `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA12, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lessthan2weeks sip_2weeks)
}


***Table A13***
clear
use "StateData.dta"
reghdfe newcases sip newcases_lag* [fw = pop], a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA13, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe newdeaths sip newdeaths_lag* [fw = pop], a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA13, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe mobility sip mobility_lag* [fw = pop], a(testingregimeid day) cluster(stateid)
qui:outreg2 using TableA13, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe unemployment sip unemployment_lag* [fw = pop], a(stateid day) cluster(stateid)
qui:outreg2 using TableA13, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)


***Table A14***
clear
use "StateData.dta"
reghdfe newdeaths sip newdeaths_lag* if p_dailydeaths != . & newdeaths_lag14 != ., a(stateid day) cluster(stateid)
qui:outreg2 using TableA14, replace dec(3)  alpha(0.01, 0.05) sortvar(sip p_dailydeaths)
reghdfe newdeaths sip p_dailydeaths if p_dailydeaths != . & newdeaths_lag14 != ., a(stateid day) cluster(stateid)
qui:outreg2 using TableA14, append dec(3)  alpha(0.01, 0.05) sortvar(sip p_dailydeaths)
reghdfe newdeaths sip p_dailydeaths newdeaths_lag* if p_dailydeaths != . & newdeaths_lag14 != ., a(stateid day) cluster(stateid)
qui:outreg2 using TableA14, append dec(3)  alpha(0.01, 0.05) sortvar(sip p_dailydeaths)
*only use the predictions from 3/26 and drop all states that had SIPs before 3/26
g sipday = day if sip == 1
egen minsip = min(sipday), by(state)
tab minsip
drop if minsip <= 22000
reghdfe newdeaths sip newdeaths_lag* if p_dailydeaths_3_26 != . & newdeaths_lag14 != ., a(stateid day) cluster(stateid)
qui:outreg2 using TableA14, append dec(3)  alpha(0.01, 0.05) sortvar(sip p_dailydeaths)
reghdfe newdeaths sip p_dailydeaths_3_26 if p_dailydeaths_3_26 != . & newdeaths_lag14 != ., a(stateid day) cluster(stateid)
qui:outreg2 using TableA14, append dec(3)  alpha(0.01, 0.05) sortvar(sip p_dailydeaths)
reghdfe newdeaths sip p_dailydeaths_3_26 newdeaths_lag* if p_dailydeaths_3_26 != . & newdeaths_lag14 != ., a(stateid day) cluster(stateid)
qui:outreg2 using TableA14, append dec(3)  alpha(0.01, 0.05) sortvar(sip p_dailydeaths)


***Table A15***
clear
use "HsiangEtAlData.dta"
reghdfe D_l_cum_confirmed_cases p_* testing_regime_*, absorb(i.adm1_id i.dow) cluster(t) //resid
outreg2 using TableA15, replace dec(3)  alpha(0.01, 0.05) 
reghdfe D_l_cum_confirmed_cases p_* /*testing_regime_change_**/, absorb(i.t i.adm1_id)  cluster(t)
outreg2 using TableA15, append dec(3)  alpha(0.01, 0.05) 
forvalues i = 1/14 {
g cases_lag`i' = D_l_cum_confirmed_cases[_n-`i'] if adm1_id == adm1_id[_n-`i'] & t == t[_n-`i'] + `i'
}
reghdfe D_l_cum_confirmed_cases p_* cases_lag1, absorb(i.t i.adm1_id)  cluster(t)
outreg2 using TableA15, append dec(3)  alpha(0.01, 0.05) 
reghdfe D_l_cum_confirmed_cases p_* cases_lag1-cases_lag7, absorb(i.t i.adm1_id)  cluster(t)
outreg2 using TableA15, append dec(3)  alpha(0.01, 0.05)


***Table A16***
clear
use "StateData.dta"
gen sip_counter = day - statesip_startday
gen sip_1_5 = 0
	replace sip_1_5 = 1 if sip_counter>=1 & sip_counter<=5
gen sip_6_14 = 0
	replace sip_6_14 = 1 if sip_counter>=6 & sip_counter<=14
gen sip_15_19 = 0
	replace sip_15_19 = 1 if sip_counter>=15 & sip_counter<=19
gen sip_20_plus = 0
	replace sip_20_plus = 1 if sip_counter>=20 & sip_counter!=.
reg logcases sip_1_5 sip_6_14 sip_15_19 sip_20_plus i.stateid##c.day i.day ///
if day>=21981 & day<=22025 & state!="NY" & state!="NJ" [fw=pop], cluster(stateid)
	outreg2 using TableA16, replace keep(sip_1_5 sip_6_14 sip_15_19 sip_20_plus) alpha(0.01, 0.05)
reg logcases sip_1_5 sip_6_14 sip_15_19 sip_20_plus i.stateid##c.day i.day ///
if day<=22025 & state!="NY" & state!="NJ" [fw=pop], cluster(stateid)
	outreg2 using TableA16, append keep(sip_1_5 sip_6_14 sip_15_19 sip_20_plus) alpha(0.01, 0.05)
reg logcases sip_1_5 sip_6_14 sip_15_19 sip_20_plus i.stateid##c.day i.day ///
if day<=22032 & state!="NY" & state!="NJ" [fw=pop], cluster(stateid)
	outreg2 using TableA16, append keep(sip_1_5 sip_6_14 sip_15_19 sip_20_plus) alpha(0.01, 0.05)
reg logcases sip_1_5 sip_6_14 sip_15_19 sip_20_plus i.stateid##c.day i.day ///
if state!="NY" & state!="NJ" [fw=pop], cluster(stateid)
	outreg2 using TableA16, append keep(sip_1_5 sip_6_14 sip_15_19 sip_20_plus) alpha(0.01, 0.05)
reg logcases sip_1_5 sip_6_14 sip_15_19 sip_20_plus i.stateid##c.day i.day [fw=pop], cluster(stateid)
	outreg2 using TableA16, append keep(sip_1_5 sip_6_14 sip_15_19 sip_20_plus) alpha(0.01, 0.05)
reg logcases sip_1_5 sip_6_14 sip_15_19 sip_20_plus i.stateid##c.day i.day, cluster(stateid)
	outreg2 using TableA16, append keep(sip_1_5 sip_6_14 sip_15_19 sip_20_plus) alpha(0.01, 0.05)
reg logcases sip_1_5 sip_6_14 sip_15_19 sip_20_plus i.stateid##c.day i.day ///
if day>=21981 & day<=22025, cluster(stateid)
	outreg2 using TableA16, append keep(sip_1_5 sip_6_14 sip_15_19 sip_20_plus) alpha(0.01, 0.05)
reg newcases sip_1_5 sip_6_14 sip_15_19 sip_20_plus i.stateid##c.day i.day ///
if day>=21981 & day<=22025 & state!="NY" & state!="NJ" [fw=pop], cluster(stateid)
	outreg2 using TableA16, append keep(sip_1_5 sip_6_14 sip_15_19 sip_20_plus) alpha(0.01, 0.05)


***Table A17***
clear
use "StateData.dta"
sort state day
forvalues i = 1/14 {
g logcases_lag`i' = logcases[_n-`i'] if state == state[_n-`i'] & day == day[_n-`i'] + `i'
}
order sip_lead7 sip_lead6 sip_lead5 sip_lead4 sip_lead3 sip_lead2 sip_lead1 sip sip_lag* logcases_lag*
foreach dv in logcases {
reghdfe `dv' sip, a(stateid day) cluster(stateid)
qui:outreg2 using TableA17, replace dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip_lead* sip , a(stateid day) cluster(stateid)
qui:outreg2 using TableA17, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip sip_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA17, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA17, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip_lead* sip sip_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA17, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip_lead* sip `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA17, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip sip_lag* `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA17, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
reghdfe `dv' sip_lead* sip sip_lag* `dv'_lag*, a(stateid day) cluster(stateid)
qui:outreg2 using TableA17, append dec(3)  alpha(0.01, 0.05) sortvar(sip_lead7-unemployment_lag7)
}



***Table A18***
clear
set more off
postutil clear
set seed 60637
use "StateData.dta"
drop stateid
egen stateid = group(state)
preserve
keep stateid day sip
save "treatments.dta", replace
restore
postfile Results outcome iteration coef using "Results.dta", replace
forvalues i = 1/1000 {
quietly {
clear
use "treatments.dta"
sort stateid
g random = uniform() if stateid != stateid[_n-1]
egen newid = group(random)
replace newid = newid[_n-1] if newid == .
drop stateid
rename newid stateid
sort stateid day
save "faketreatments.dta", replace
clear
use "StateData.dta"
drop sip
sort stateid day
merge stateid day using "faketreatments.dta"
tab _merge
drop _merge
reghdfe newcases sip newcases_lag*, a(stateid day) cluster(stateid)
post Results (1) (`i') (_b[sip]) 
reghdfe newdeaths sip newdeaths_lag*, a(stateid day) cluster(stateid)
post Results (2) (`i') (_b[sip]) 
reghdfe mobility sip mobility_lag*, a(stateid day) cluster(stateid)
post Results (3) (`i') (_b[sip]) 
reghdfe unemployment sip unemployment_lag*, a(stateid day) cluster(stateid)
post Results (4) (`i') (_b[sip]) 
}
disp `i'
}
postclose Results
clear
set more off
use "StateData.dta"
reghdfe newcases sip newcases_lag*, a(stateid day) cluster(stateid)
local coef1 = _b[sip]
reghdfe newdeaths sip newdeaths_lag*, a(stateid day) cluster(stateid)
local coef2 = _b[sip]
reghdfe mobility sip mobility_lag*, a(stateid day) cluster(stateid)
local coef3 = _b[sip]
reghdfe unemployment sip unemployment_lag*, a(stateid day) cluster(stateid)
local coef4 = _b[sip]
clear
use "Results.dta"
forvalues i = 1/4 {
preserve
keep if outcome == `i'
sum coef
sort coef
sum coef if _n == 26
sum coef if _n == 975
g greater = abs(coef) > abs(`coef`i'')
sum greater
restore
}


***Figure A1***
clear
set more off
use "StateData.dta"
reghdfe newcases sip newcases_lag*, a(stateid day) cluster(stateid)
local coef1 = _b[sip]
reghdfe newdeaths sip newdeaths_lag*, a(stateid day) cluster(stateid)
local coef2 = _b[sip]
reghdfe mobility sip mobility_lag*, a(stateid day) cluster(stateid)
local coef3 = _b[sip]
reghdfe unemployment sip unemployment_lag*, a(stateid day) cluster(stateid)
local coef4 = _b[sip]
clear
use "Results.dta"
forvalues i = 1/4 {
preserve
keep if outcome == `i'
sum coef
sort coef
sum coef if _n == 26
sum coef if _n == 975
g greater = abs(coef) > abs(`coef`i'')
sum greater
restore
egen maxsip = max(sip), by(state)
g never = maxsip < 1
g sipday = day if sip == 1
egen maxsipday = max(sipday), by(state)
g stayon = maxsipday == 22066
g onoff = never + stayon == 0
drop max*
*recode day as Days after Feb 24
replace day = day - 21969
foreach i in AL AR AZ CA CO CT DC DE FL GA HI IA ID IL IN KS KY LA MA MD ME MI MN ///
MO MS MT NC ND NE NH NJ NM NV NY OH OK OR PA RI SC SD TN TX UT VA VT WA WI WV WY {
preserve
keep if state == "`i'"
replace newcases = newcases/100
lpoly newcases day, bw(7) nograph gen(k_newcases) at(day)
replace newcases = . if newcases > 5
sum never
if r(mean) == 1 {
	sort day
	graph twoway (line k_newcases day) (scatter newcases day)
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush New Cases
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 5 1 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore
	}
sum stayon
if r(mean) == 1 {	
	sum day if sip == 1
	local start = r(min)
	sort day
	graph twoway (line k_newcases day) (scatter newcases day) , xline(`start')
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush New Cases
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 5 1 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore	
	}
sum onoff
if r(mean) == 1 {
	sum day if sip == 1
	local start = r(min)
	local end = r(max)
	sort day
	graph twoway (line k_newcases day) (scatter newcases day) , xline(`start') xline(`end')
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush New Cases
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 5 1 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
	gr_edit .plotregion1._xylines[2].style.editstyle linestyle(color(black) pattern(dot)) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore
	}	
}
graph combine "AL.gph" "AR.gph" "AZ.gph" "CA.gph" "CO.gph" "CT.gph" "DC.gph" ///
"DE.gph" "FL.gph" "GA.gph" "HI.gph" "IA.gph" "ID.gph" "IL.gph" "IN.gph" "KS.gph" ///
"KY.gph" "LA.gph" "MA.gph" "MD.gph" "ME.gph" "MI.gph" "MN.gph" "MO.gph" ///
"MS.gph" "MT.gph" "NC.gph" "ND.gph" "NE.gph" "NH.gph" "NJ.gph" "NM.gph" ///
"NV.gph" "NY.gph" "OH.gph" "OK.gph" "OR.gph" "PA.gph" "RI.gph" "SC.gph" ///
"SD.gph" "TN.gph" "TX.gph" "UT.gph" "VA.gph" "VT.gph" "WA.gph" "WI.gph" ///
"WV.gph" "WY.gph", col(5)
gr_edit .style.editstyle declared_ysize(7) editcopy
gr_edit .style.editstyle margin(zero) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
graph export "FigureA1.png", replace as(png)


***Figure A2***
clear
set more off
use "StateData.dta"
egen maxsip = max(sip), by(state)
g never = maxsip < 1
g sipday = day if sip == 1
egen maxsipday = max(sipday), by(state)
g stayon = maxsipday == 22066
g onoff = never + stayon == 0
drop max*
*recode day as Days after Feb 24
replace day = day - 21969
foreach i in AL AR AZ CA CO CT DC DE FL GA HI IA ID IL IN KS KY LA MA MD ME MI MN ///
MO MS MT NC ND NE NH NJ NM NV NY OH OK OR PA RI SC SD TN TX UT VA VT WA WI WV WY {
preserve
keep if state == "`i'"
replace newdeaths = newdeaths/10
lpoly newdeaths day, bw(7) nograph gen(k_newdeaths) at(day)
replace newdeaths = . if newdeaths > 6
sum never
if r(mean) == 1 {
	sort day
	graph twoway (line k_newdeaths day) (scatter newdeaths day)
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush New Cases
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 6 2 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore
	}
sum stayon
if r(mean) == 1 {	
	sum day if sip == 1
	local start = r(min)
	sort day
	graph twoway (line k_newdeaths day) (scatter newdeaths day) , xline(`start')
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush New Cases
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 6 2 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore	
	}
sum onoff
if r(mean) == 1 {
	sum day if sip == 1
	local start = r(min)
	local end = r(max)
	sort day
	graph twoway (line k_newdeaths day) (scatter newdeaths day) , xline(`start') xline(`end')
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush Deaths
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 6 2 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
	gr_edit .plotregion1._xylines[2].style.editstyle linestyle(color(black) pattern(dot)) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore
	}	
}
graph combine "AL.gph" "AR.gph" "AZ.gph" "CA.gph" "CO.gph" "CT.gph" "DC.gph" ///
"DE.gph" "FL.gph" "GA.gph" "HI.gph" "IA.gph" "ID.gph" "IL.gph" "IN.gph" "KS.gph" ///
"KY.gph" "LA.gph" "MA.gph" "MD.gph" "ME.gph" "MI.gph" "MN.gph" "MO.gph" ///
"MS.gph" "MT.gph" "NC.gph" "ND.gph" "NE.gph" "NH.gph" "NJ.gph" "NM.gph" ///
"NV.gph" "NY.gph" "OH.gph" "OK.gph" "OR.gph" "PA.gph" "RI.gph" "SC.gph" ///
"SD.gph" "TN.gph" "TX.gph" "UT.gph" "VA.gph" "VT.gph" "WA.gph" "WI.gph" ///
"WV.gph" "WY.gph", col(5)
gr_edit .style.editstyle declared_ysize(7) editcopy
gr_edit .style.editstyle margin(zero) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
graph export "FigureA2.png", replace as(png)


***Figure A3***
clear
set more off
use "StateData.dta"
egen maxsip = max(sip), by(state)
g never = maxsip < 1
g sipday = day if sip == 1
egen maxsipday = max(sipday), by(state)
g stayon = maxsipday == 22066
g onoff = never + stayon == 0
drop max*
*recode day as Days after Feb 24
replace day = day - 21969
foreach i in AL AR AZ CA CO CT DC DE FL GA HI IA ID IL IN KS KY LA MA MD ME MI MN ///
MO MS MT NC ND NE NH NJ NM NV NY OH OK OR PA RI SC SD TN TX UT VA VT WA WI WV WY {
preserve
keep if state == "`i'"
lpoly unemployment day, bw(14) nograph gen(k_unemployment) at(day)
sum never
if r(mean) == 1 {
	sort day
	graph twoway (line k_unemployment day) (scatter unemployment day)
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush Unemployment
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 30 10 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore
	}
sum stayon
if r(mean) == 1 {	
	sum day if sip == 1
	local start = r(min)
	sort day
	graph twoway (line k_unemployment day) (scatter unemployment day) , xline(`start')
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush Unemployment
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 30 10 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore	
	}
sum onoff
if r(mean) == 1 {
	sum day if sip == 1
	local start = r(min)
	local end = r(max)
	sort day
	graph twoway (line k_unemployment day) (scatter unemployment day) , xline(`start') xline(`end')
	gr_edit .legend.draw_view.setstyle, style(no)
	gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
	gr_edit .yaxis1.title.text = {}
	gr_edit .yaxis1.title.text.Arrpush Unemployment
	gr_edit .xaxis1.title.text = {}
	gr_edit .xaxis1.title.text.Arrpush Day
	gr_edit .xaxis1.reset_rule 0 100 20 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.reset_rule 0 30 10 , tickset(major) ruletype(range) 
	gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
	gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
	gr_edit .plotregion1._xylines[2].style.editstyle linestyle(color(black) pattern(dot)) editcopy
	gr_edit .plotregion1.plot1.style.editstyle line(color(gs8) width(thick)) editcopy
	gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(vsmall)) editcopy
	gr_edit .title.style.editstyle color(black) editcopy
	gr_edit .title.text = {}
	gr_edit .title.text.Arrpush `i'
	graph save "`i'.gph", replace
	restore
	}	
}
graph combine "AL.gph" "AR.gph" "AZ.gph" "CA.gph" "CO.gph" "CT.gph" "DC.gph" ///
"DE.gph" "FL.gph" "GA.gph" "HI.gph" "IA.gph" "ID.gph" "IL.gph" "IN.gph" "KS.gph" ///
"KY.gph" "LA.gph" "MA.gph" "MD.gph" "ME.gph" "MI.gph" "MN.gph" "MO.gph" ///
"MS.gph" "MT.gph" "NC.gph" "ND.gph" "NE.gph" "NH.gph" "NJ.gph" "NM.gph" ///
"NV.gph" "NY.gph" "OH.gph" "OK.gph" "OR.gph" "PA.gph" "RI.gph" "SC.gph" ///
"SD.gph" "TN.gph" "TX.gph" "UT.gph" "VA.gph" "VT.gph" "WA.gph" "WI.gph" ///
"WV.gph" "WY.gph", col(5)
gr_edit .style.editstyle declared_ysize(7) editcopy
gr_edit .style.editstyle margin(zero) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
graph export "FigureA3.png", replace as(png)


***Figure A4***
clear
use "StateData.dta"
sort state day
forvalues i = 1/14 {
g logcases_lag`i' = logcases[_n-`i'] if state == state[_n-`i'] & day == day[_n-`i'] + `i'
}
reghdfe logcases sip sip_lag* logcases_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Log Cases
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 14 2 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph export "FigureA4.png", replace as(png)


***Figure A5***
clear
postutil clear
set more off
use "StateData.dta"
drop *lag* *lead*
g previousdays = 0
sort state day
forvalues i = 1/21 {
g sip_lag`i' = sip[_n-`i'] if state == state[_n-`i'] & day == day[_n-`i'] + `i'
replace previousdays = previousdays + 1 if sip_lag`i' == 1
}
forvalues i = 1/7 {
g sip_lead`i' = sip[_n+`i'] if state == state[_n+`i'] & day == day[_n+`i'] - `i'
}
forvalues i = 1/14 {
g newcases_lag`i' = newcases[_n-`i'] if state == state[_n-`i'] & day == day[_n-`i'] + `i'
g newdeaths_lag`i' = newdeaths[_n-`i'] if state == state[_n-`i'] & day == day[_n-`i'] + `i'
}
forvalues i = 1/14 {
g mobility_lag`i' = mobility[_n-`i'] if state == state[_n-`i'] & day == day[_n-`i'] + `i'
}
g unemployment_lag7 = unemployment[_n-7] if state == state[_n-7] & day == day[_n-7] + 7
g unemployment_lag14 = unemployment[_n-14] if state == state[_n-14] & day == day[_n-14] + 14
save "StateLevelData_21.dta", replace
*Newcases
clear
use "StateLevelData_21.dta"
reghdfe newcases sip sip_lag* newcases_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
scalar effect = effect + _b[sip_lag15]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 = 0
post Results (15) (effect) (r(p))
scalar effect = effect + _b[sip_lag16]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 = 0
post Results (16) (effect) (r(p))
scalar effect = effect + _b[sip_lag17]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 = 0
post Results (17) (effect) (r(p))
scalar effect = effect + _b[sip_lag18]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 = 0
post Results (18) (effect) (r(p))
scalar effect = effect + _b[sip_lag19]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 = 0
post Results (19) (effect) (r(p))
scalar effect = effect + _b[sip_lag20]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 + sip_lag20 = 0
post Results (20) (effect) (r(p))
scalar effect = effect + _b[sip_lag21]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 + sip_lag20 + sip_lag21  = 0
post Results (21) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on New Cases per Million Residents
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 21 3 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 21 `"21+"', tickset(major)
graph save "LaggedEffects_NewCases.gph", replace
*Deaths
clear
use "StateLevelData_21.dta"
reghdfe newdeaths sip sip_lag* newdeaths_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
scalar effect = effect + _b[sip_lag15]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 = 0
post Results (15) (effect) (r(p))
scalar effect = effect + _b[sip_lag16]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 = 0
post Results (16) (effect) (r(p))
scalar effect = effect + _b[sip_lag17]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 = 0
post Results (17) (effect) (r(p))
scalar effect = effect + _b[sip_lag18]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 = 0
post Results (18) (effect) (r(p))
scalar effect = effect + _b[sip_lag19]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 = 0
post Results (19) (effect) (r(p))
scalar effect = effect + _b[sip_lag20]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 + sip_lag20 = 0
post Results (20) (effect) (r(p))
scalar effect = effect + _b[sip_lag21]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 + sip_lag20 + sip_lag21  = 0
post Results (21) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Deaths per Million Residents
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 21 3 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 21 `"21+"', tickset(major)
graph save "LaggedEffects_Deaths.gph", replace
*Mobility
clear
use "StateLevelData_21.dta"
reghdfe mobility sip sip_lag* mobility_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
scalar effect = effect + _b[sip_lag15]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 = 0
post Results (15) (effect) (r(p))
scalar effect = effect + _b[sip_lag16]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 = 0
post Results (16) (effect) (r(p))
scalar effect = effect + _b[sip_lag17]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 = 0
post Results (17) (effect) (r(p))
scalar effect = effect + _b[sip_lag18]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 = 0
post Results (18) (effect) (r(p))
scalar effect = effect + _b[sip_lag19]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 = 0
post Results (19) (effect) (r(p))
scalar effect = effect + _b[sip_lag20]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 + sip_lag20 = 0
post Results (20) (effect) (r(p))
scalar effect = effect + _b[sip_lag21]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 + sip_lag20 + sip_lag21  = 0
post Results (21) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Mobility
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 21 3 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 21 `"21+"', tickset(major)
graph save "LaggedEffects_Mobility.gph", replace
*Unemployment
clear
use "StateLevelData_21.dta"
reghdfe unemployment sip sip_lag* unemployment_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip]
test sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
scalar effect = effect + _b[sip_lag15]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 = 0
post Results (15) (effect) (r(p))
scalar effect = effect + _b[sip_lag16]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 = 0
post Results (16) (effect) (r(p))
scalar effect = effect + _b[sip_lag17]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 = 0
post Results (17) (effect) (r(p))
scalar effect = effect + _b[sip_lag18]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 = 0
post Results (18) (effect) (r(p))
scalar effect = effect + _b[sip_lag19]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 = 0
post Results (19) (effect) (r(p))
scalar effect = effect + _b[sip_lag20]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 + sip_lag20 = 0
post Results (20) (effect) (r(p))
scalar effect = effect + _b[sip_lag21]
test sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14 + sip_lag15 + sip_lag16 + sip_lag17 + sip_lag18 + sip_lag19 + sip_lag20 + sip_lag21  = 0
post Results (21) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Insured Unemployment Rate
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule 0 21 3 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 21 `"21+"', tickset(major)
graph save "LaggedEffects_Unemployment.gph", replace
*Combining
graph combine "LaggedEffects_NewCases.gph" "LaggedEffects_Deaths.gph" "LaggedEffects_Mobility.gph" "LaggedEffects_Unemployment.gph", col(2)
gr_edit .plotregion1.graph1.yaxis1.title.text = {}
gr_edit .plotregion1.graph1.yaxis1.title.text.Arrpush Effect of SIP
gr_edit .plotregion1.graph3.yaxis1.title.text = {}
gr_edit .plotregion1.graph3.yaxis1.title.text.Arrpush Effect of SIP
gr_edit .plotregion1.graph1.xaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph2.xaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph2.yaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph4.yaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph1.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph1.title.text = {}
gr_edit .plotregion1.graph1.title.text.Arrpush Cases
gr_edit .plotregion1.graph2.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph2.title.text = {}
gr_edit .plotregion1.graph2.title.text.Arrpush Deaths
gr_edit .plotregion1.graph3.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph3.title.text = {}
gr_edit .plotregion1.graph3.title.text.Arrpush Mobility
gr_edit .plotregion1.graph4.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph4.title.text = {}
gr_edit .plotregion1.graph4.title.text.Arrpush Unemployment
gr_edit .style.editstyle margin(zero) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
graph export "FigureA5.png", replace as(png)


***Figure A6***
*Newcases
clear
use "StateData.dta"
reghdfe newcases sip_lead* sip sip_lag* newcases_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip_lead7]
test sip_lead7 = 0
post Results (-7) (effect) (r(p))
scalar effect = _b[sip_lead6]
test sip_lead7 + sip_lead6 = 0
post Results (-6) (effect) (r(p))
scalar effect = _b[sip_lead5]
test sip_lead7 + sip_lead6 + sip_lead5 = 0
post Results (-5) (effect) (r(p))
scalar effect = _b[sip_lead4]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 = 0
post Results (-4) (effect) (r(p))
scalar effect = _b[sip_lead3]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 = 0
post Results (-3) (effect) (r(p))
scalar effect = _b[sip_lead2]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 = 0
post Results (-2) (effect) (r(p))
scalar effect = _b[sip_lead1]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 = 0
post Results (-1) (effect) (r(p))
scalar effect = _b[sip]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0) xline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on New Cases per Million Residents
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule -7 14 3, tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1._xylines[2].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph save "LaggedEffects_NewCases.gph", replace
*Deaths
clear
use "StateData.dta"
reghdfe newdeaths sip_lead* sip sip_lag* newdeaths_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip_lead7]
test sip_lead7 = 0
post Results (-7) (effect) (r(p))
scalar effect = _b[sip_lead6]
test sip_lead7 + sip_lead6 = 0
post Results (-6) (effect) (r(p))
scalar effect = _b[sip_lead5]
test sip_lead7 + sip_lead6 + sip_lead5 = 0
post Results (-5) (effect) (r(p))
scalar effect = _b[sip_lead4]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 = 0
post Results (-4) (effect) (r(p))
scalar effect = _b[sip_lead3]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 = 0
post Results (-3) (effect) (r(p))
scalar effect = _b[sip_lead2]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 = 0
post Results (-2) (effect) (r(p))
scalar effect = _b[sip_lead1]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 = 0
post Results (-1) (effect) (r(p))
scalar effect = _b[sip]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0) xline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Deaths per Million Residents
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule -7 14 3 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1._xylines[2].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph save "LaggedEffects_Deaths.gph", replace
*Mobility
clear
use "StateData.dta"
reghdfe mobility sip_lead* sip sip_lag* mobility_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip_lead7]
test sip_lead7 = 0
post Results (-7) (effect) (r(p))
scalar effect = _b[sip_lead6]
test sip_lead7 + sip_lead6 = 0
post Results (-6) (effect) (r(p))
scalar effect = _b[sip_lead5]
test sip_lead7 + sip_lead6 + sip_lead5 = 0
post Results (-5) (effect) (r(p))
scalar effect = _b[sip_lead4]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 = 0
post Results (-4) (effect) (r(p))
scalar effect = _b[sip_lead3]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 = 0
post Results (-3) (effect) (r(p))
scalar effect = _b[sip_lead2]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 = 0
post Results (-2) (effect) (r(p))
scalar effect = _b[sip_lead1]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 = 0
post Results (-1) (effect) (r(p))
scalar effect = _b[sip]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0) xline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Mobility
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule -7 14 3 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1._xylines[2].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph save "LaggedEffects_Mobility.gph", replace
*Unemployment
clear
use "StateData.dta"
reghdfe unemployment sip_lead* sip sip_lag* unemployment_lag*, a(stateid day) cluster(stateid)
postfile Results days effect pval using "Results.dta", replace
scalar effect = _b[sip_lead7]
test sip_lead7 = 0
post Results (-7) (effect) (r(p))
scalar effect = _b[sip_lead6]
test sip_lead7 + sip_lead6 = 0
post Results (-6) (effect) (r(p))
scalar effect = _b[sip_lead5]
test sip_lead7 + sip_lead6 + sip_lead5 = 0
post Results (-5) (effect) (r(p))
scalar effect = _b[sip_lead4]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 = 0
post Results (-4) (effect) (r(p))
scalar effect = _b[sip_lead3]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 = 0
post Results (-3) (effect) (r(p))
scalar effect = _b[sip_lead2]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 = 0
post Results (-2) (effect) (r(p))
scalar effect = _b[sip_lead1]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 = 0
post Results (-1) (effect) (r(p))
scalar effect = _b[sip]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip = 0
post Results (0) (effect) (r(p))
scalar effect = effect + _b[sip_lag1]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 = 0
post Results (1) (effect) (r(p))
scalar effect = effect + _b[sip_lag2]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 = 0
post Results (2) (effect) (r(p))
scalar effect = effect + _b[sip_lag3]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 = 0
post Results (3) (effect) (r(p))
scalar effect = effect + _b[sip_lag4]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 = 0
post Results (4) (effect) (r(p))
scalar effect = effect + _b[sip_lag5]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 = 0
post Results (5) (effect) (r(p))
scalar effect = effect + _b[sip_lag6]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 = 0
post Results (6) (effect) (r(p))
scalar effect = effect + _b[sip_lag7]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 = 0
post Results (7) (effect) (r(p))
scalar effect = effect + _b[sip_lag8]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 = 0
post Results (8) (effect) (r(p))
scalar effect = effect + _b[sip_lag9]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 = 0
post Results (9) (effect) (r(p))
scalar effect = effect + _b[sip_lag10]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 = 0
post Results (10) (effect) (r(p))
scalar effect = effect + _b[sip_lag11]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 = 0
post Results (11) (effect) (r(p))
scalar effect = effect + _b[sip_lag12]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12  = 0
post Results (12) (effect) (r(p))
scalar effect = effect + _b[sip_lag13]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 = 0
post Results (13) (effect) (r(p))
scalar effect = effect + _b[sip_lag14]
test sip_lead7 + sip_lead6 + sip_lead5 + sip_lead4 + sip_lead3 + sip_lead2 + sip_lead1 + sip + sip_lag1 + sip_lag2 + sip_lag3 + sip_lag4 + sip_lag5 + sip_lag6 + sip_lag7 + sip_lag8 + sip_lag9 + sip_lag10 + sip_lag11 + sip_lag12 + sip_lag13 + sip_lag14  = 0
post Results (14) (effect) (r(p))
postclose Results
clear
use "Results.dta"
g effect_significant = effect if pval < .05
g effect_insignificant = effect if pval >= .05
sort days
graph twoway (line effect days) (scatter effect_* days), yline(0) xline(0)
gr_edit .legend.draw_view.setstyle, style(no)
gr_edit .style.editstyle margin(vsmall) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
gr_edit .yaxis1.title.text = {}
gr_edit .yaxis1.title.text.Arrpush Effect of SIP on Insured Unemployment Rate
gr_edit .xaxis1.title.text = {}
gr_edit .xaxis1.title.text.Arrpush Days SIP Has Been in Place
gr_edit .xaxis1.reset_rule -7 14 3 , tickset(major) ruletype(range) 
gr_edit .yaxis1.style.editstyle majorstyle(gridstyle(linestyle(color(white)))) editcopy
gr_edit .plotregion1._xylines[1].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1._xylines[2].style.editstyle linestyle(color(black) pattern(dash)) editcopy
gr_edit .plotregion1.plot1.style.editstyle line(color(black) width(thick)) editcopy
gr_edit .plotregion1.plot3.style.editstyle marker(fillcolor(black) linestyle(color(black)) size(medlarge)) editcopy
gr_edit .plotregion1.plot2.style.editstyle marker(fillcolor(red) linestyle(color(red)) size(medlarge)) editcopy
gr_edit .xaxis1.major.num_rule_ticks = 8
gr_edit .xaxis1.edit_tick 8 14 `"14+"', tickset(major)
graph save "LaggedEffects_Unemployment.gph", replace
graph combine "LaggedEffects_NewCases.gph" "LaggedEffects_Deaths.gph" "LaggedEffects_Mobility.gph" "LaggedEffects_Unemployment.gph", col(2)
gr_edit .plotregion1.graph1.yaxis1.title.text = {}
gr_edit .plotregion1.graph1.yaxis1.title.text.Arrpush Effect of SIP
gr_edit .plotregion1.graph3.yaxis1.title.text = {}
gr_edit .plotregion1.graph3.yaxis1.title.text.Arrpush Effect of SIP
gr_edit .plotregion1.graph1.xaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph2.xaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph2.yaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph4.yaxis1.title.draw_view.setstyle, style(no)
gr_edit .plotregion1.graph1.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph1.title.text = {}
gr_edit .plotregion1.graph1.title.text.Arrpush Cases
gr_edit .plotregion1.graph2.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph2.title.text = {}
gr_edit .plotregion1.graph2.title.text.Arrpush Deaths
gr_edit .plotregion1.graph3.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph3.title.text = {}
gr_edit .plotregion1.graph3.title.text.Arrpush Mobility
gr_edit .plotregion1.graph4.title.style.editstyle color(black) editcopy
gr_edit .plotregion1.graph4.title.text = {}
gr_edit .plotregion1.graph4.title.text.Arrpush Unemployment
gr_edit .style.editstyle margin(zero) boxstyle(shadestyle(color(white)) linestyle(color(white))) editcopy
graph export "TableA6.png", replace as(png)
